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ABSTRACT 

We explore the presence of torques at the inner edges of geometrically-thin black hole accretion disks using 
3 -dimensional magnetohydrodynamic (MHD) simulations in a pseudo-Newtonian potential. By varying the sat- 
uration level of the magnetorotational instability that leads to angular momentum transport, we show that the 
dynamics of gas inside the radius of marginal stability varies depending upon the magnetic field strength just out- 
side that radius. Weak fields are unable to causally connect material within the plunging region to the rest of the 
disk, and zero torque is an approximately correct boundary condition at the radius of marginal stability. Stronger 
fields, which we obtain artificially but which may occur physically within more complete disk models, are able 
to couple at least some parts of the plunging region to the rest of the disk. In this case, angular momentum (and 
implicitly energy) is extracted from the material in the plunging region. Furthermore, the magnetic coupling to 
the plunging region can be highly time dependent with large fluctuations in the torque at the radius of marginal 
stability. This implies varying accretion efficiencies, both across systems and within a given system at different 
times. The results suggest a possible link between changes in X-ray and outflow activity, with both being driven 
by transitions between weak and strong field states. 

Subject headings: accretion, accretion disks - black hole physics - MHD - hydrodynamics - instabilities 



1. INTRODUCTION 

In a wide class of black hole systems, including QSOs, 
Seyfert galaxies, and high-state Galactic Black Hole Candi- 
, dates, the accretion disk is almost certainly geometrically thin 
■ and radiatively efficient. The radiation emitted from an annu- 
lus of such a disk does not simply equal the change in grav- 
. itational binding energy as gas flows across the annulus, but 
' rather includes a contribution from energy transported into the 
annulus from elsewhere in the disk. As a consequence, the disk 
structure, along with related quantities such as the radiative effi- 
ciency (the fraction of rest mass energy which is radiated during 
the accretion process), depends upon the location and nature of 
the boundary conditions at the inner disk edge. 

For a black hole accretion disk, a natural location to place 
the inner boundary condition is at the radius of marginal sta- 
bility, r = Tnis. This is the radius inside of which circu- 
lar orbits are no longer stable, and is at 6GM/c^ for a non- 
rotating (Schwarzschild) black hole. In the standard model for 
black hole accretion (Novikov & Thorne 1973; Paczynski & 
Bisnovatyi-Kogan 1981; Abramowicz & Kato 1989; Paczynski 
2000), one assumes that there is no angular momentum trans- 
port across r = r^ns- The motivation for this zero torque bound- 
ary condition (ZTBC) is that the radial flow rapidly becomes su- 
personic once inside r — rms, and so the material loses causal 
contact with the rest of the disk. Thereafter, material simply 
spirals ballistically into the black hole. For this reason, we will 
refer to the region r < as the plunging region. 

The ZTBC limits the accretion efficiency to around 6% (for 
a non-rotating black hole) - already a substantial figure. How- 
ever, it is possible that even higher efficiencies may arise. The 



suggestion is that magnetic fields, which are generated in the 
disk by the magnetorotational instability (Balbus & Hawley 
1991) may become dynamically significant within the plung- 
ing region. These fields could then causally connect the plung- 
ing region to the rest of the accretion disk, allowing energy and 
angular momentum to be extracted from gas as it executes its 
final spiral into the black hole (Krolik 1999; Gammie 1999; 
Agol & Krolik 2000). Subsequent numerical simulations (none 
of which are General Relativistic) have confirmed the presence 
of magnetic torques at the marginally stable orbit, but have 
reached different conclusions as to their impact on the dynam- 
ics of the flow. Global MHD simulations of adiabatic accretion 
tori (Hawley 2000; Hawley & Krolik 2001) showed outward 
angular momentum transport continuing within r^ns, while un- 
stratified simulations suggested more modest effects (Armitage, 
Reynolds & Chiang 2001; Hawley 2001). The important pa- 
rameters or numerical differences that cause this different be- 
havior remain to be determined. 

In this letter we explore the suggestion (Charles Gammie, 
private communication) that the dynamics of the flow within 
the plunging region may depend upon the strength of magnetic 
fields in the disk at r > r,ns. We extend our earlier work by 
presenting 'boosted' simulations in which the strength of the 
magnetic field in the saturated turbulence is enhanced by start- 
ing with an initially large seed field. This is a numerical trick, 
though it may have a physical counterpart if the inner disk is 
threaded by a large-scale magnetic field. For the purposes of 
this letter, however, the aim is solely to allow us to study how 
the dynamics of the plunging region vary with field strength 
in a controlled fashion. We show that the dynamics of gas in 
the plunging region, and the validity of the ZTBC, do indeed 
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depend upon the saturation field strength, and discuss the im- 
plications for observations of accreting black holes. 

2. THE SIMULATIONS 

We use the ZEUS code to solve the equations of ideal MHD 
(Stone & Norman 1992a, 1992b; Clarke, Norman & Fielder 
1994; Norman 2000) in cylindrical co-ordinates, using a setup 
very similar to that described previously (Armitage, Reynolds 
& Chiang 2001). Stated briefly, the computational domain 
is flie wedge given by r G (ri„,ro,ii), (t> G {0,(f>^i,^) and 
z G (— -Zmax, -Zmax)- The cffccts of General Relativity (in 
particular the existence of an innermost stable orbit) are mod- 
elled in this inherently non-relativistic code by using a pseudo- 
Newtonian potential (Paczynski & Wiita 1980), $ = - ^^i^ , 

' ' Sch 

where rsch = 2GM/(?. The vertical component of gravity is 
neglected, so the disk has no vertical structure. There is there- 
fore no time-averaged variation of density or magnetic field 
strength with z. The equation of state is isothermal, with a 
sound speed Cg which is assumed uniform and constant. At 
rms> the ratio of the sound speed to the Keplerian velocity is 
Cs/i^t^ = 0.065, slightly cooler than our previous simulations. 

The initial conditions for density and velocity were obtained 
by relaxing a Gaussian radial density profile, using ZEUS in 
its ID hydrodynamic mode, to produce a density profile which 
is in accurate numerical equilibrium. The initial configurations 
of the 3D MHD simulations were then produced by adding a 
seed magnetic field. Table 1 shows the computational domain, 
resolution, initial value of /3, and seed field configuration for 
the simulations. As usual, (3 is defined as the ratio of thermal 
energy density to magnetic energy density. In all cases, we set 
r,„ = 4GAf/c2 and ro„i = 20GM/c^. The radial bound- 
ary conditions were such as to permit outflow, and all other 
boundaries were made periodic. These simulations were then 
evolved until the magnetorotational instability (MRl) became 
non-linear and produced fully developed MHD turbulence (in 
the sense that the MRI has saturated and there is no further long- 
term growth of the magnetic energy). This turbulence drives 
accretion of material through the inner radial boundary, and the 
run is continued until an appreciable amount of material has left 
the computational domain.. 

3. THE DYNAMICS OF THE FLOW WITHIN THE PLUNGING REGION 

3.1. Time averaged properties 

We consider first the time-averaged properties of our runs. 
Figure 1 shows time-averaged radial profiles for (i.e., the 
ratio of magnetic field energy density to thermal energy density) 
and the specific angular momentum I. As intended, the mag- 
netic field in Run 2, which starts with a relatively strong ver- 
tical seed field (/J^ = 500), saturates with significantly higher 
relative magnetic field energy than those runs which start with 
either a weaker initial field or a comparable toroidal field. This 
is due to the continued boosting of the MRI from the conserved 
vertical flux that accompanies the initial vertical field (Hawley, 
Gammie & Balbus 1995). Increasing the size of the computa- 
tional domain also increases the saturation value of the mag- 
netic field, but by a smaUer factor. 

The dynamics of the flow - diagnosed using the specific an- 
gular momentum - correlates with the saturation field strength. 
Firstly, consider the unboosled runs (Runs 1,3 and 4; shown by 
the solid, dashed and dot-dashed lines respectively). As one 
approaches and enters the plunging region, the radial profile of 
I flattens. If the ZTBC is strictly vahd, I should be constant 



within the plunging region. For these runs, we obtain the same 
result as before (Armitage, Reynolds & Chiang 2001), namely 
that the ZTBC yields a good but not perfect description of the 
inner disk dynamics. Indeed, examination of the r(j) component 
of the magnetic stress tensor shows that the Maxwell stresses at 
T = ?'ms do not vanish but are less than the peak stress, which 
occurs at r 10 - 12 GMjc?, by a factor of 5-8. 

On the other hand, angular momentum transport within the 
plunging region is much more significant in the boosted run 
(Run 2; dotted line in Fig. 1). The time averaged stress at 
T = Ttas is Still less than the peak value (which occurs at 
r ~ 8 GMj c^), but only by a factor of 2. Significant stresses 
continue to operate within the region r < rms, and there is evi- 
dently a clear violation of the ZTBC. 

3.2. Variability within individual simulations 

Turbulent accretion disks are inherently variable systems and 
so it is also interesting to examine variability during our simula- 
tions. For this purpose, we will compare and contrast Run 1 and 
Run 2. Both runs display order unity temporal and spatial fluc- 
tuations in For the unboosted run these fluctuations have 
only a small effect on /. For the boosted run shown in Fig. 2, 
however, there are much larger fluctuations in L In particularly, 
during one timeslice from Run 2, a large part of the plunging 
region becomes magnetically dominated (/3 < 1). During this 
period (shown with a thick line in the bottom panels of Fig. 2) 
I declines within the plunging region with the same slope as 
in the rest of the disk, with the Maxwell stresses staying fairly 
constant well into the plunging region. Events as dramatic as 
this are fairly rare (caught in 1 out of the 12 independent times- 
lices that we analyzed), but it is clear that the stress at r = 
is a rapidly varying function of time, with much more frequent 
occurrences of moderate stress. 

Extraction of energy and angular momentum from within 
the plunging region can only occur if the material can remain 
causally connected to the rest of the disk. Analysis of the sim- 
ulations shows that large stresses occur within the plunging re- 
gion when the magnetic fields are strong enough for the peak 
radial MHD wave speed to exceed the radial inflow velocity 
throughout a large part of the plunging region. It is worth not- 
ing that even during the magnetically dominated event seen in 
Run 2, the azimuthally averaged MHD wave speeds are always 
slower that the inflow speeds. The angular momentum transport 
within this highly inhomogeneous plasma therefore appears to 
be mediated by low-/3 (i.e., higher relative magnetic field) fila- 
ments. 

4. DISCUSSION 

4. 1 . Radiative efficiency 

Calculating the implied change to the radiative efficiency as 
a consequence of the magnetic torque at r^s is not straightfor- 
ward (see the discussion in Hawley & Krolik 2001). A simple 
approach is to note that the dissipation rate is 

where T'"'^ is the local stress. We can then estimate the change 
in the radiative efficiency by comparing Q{r), computed using 
the stress measured in the simulations, with the form expected 
in a steady disk if there were no stress at rms- 

In practice, the non-steady nature of the flow, and the lim- 
ited radial extent of the simulations, creates difficulties. An 
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Table 1 

The set of simulations. Columns show (1) Run number, (2) 0-domain, (3) ^-domain, (4) number of radial cells (and 
number inside r = tms), (5) number of (j) cells, (6) number of vertical cells, and (7) initial /3 of plasma (and direction). 
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Fig. 1 . — Time-averaged radial profiles of /3 ^ (the ratio of magnetic energy density to thermal energy density) and I (specific angular momentum) for our high 
resolution simulations. The runs are denoted by: Run 1 (solid line), Run 2 (dotted line). Run 3 (dashed line), Run 4 (dot-dashed line). 
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Fig. 2. — Radial profiles of /3 ^ and I for individual timeslices from the boosted simulation. Run 2. In particular, note the anomalous event (shown with a thick 
Une) during which the plimging region becomes magnetically dominated and there is efficient angular momentum transport out of this region. 



indicative estimate, however, follows from noting that in run 2 
the mean magnetic stress is roughly constant between and 
2rms- Compared to a standard disc model, this implies a 50% 
increase in the dissipation in the region < r < 2rms> and 
a 10% increase over < r < lOrms- Significantiy larger 
effects are seen in some individual timeslices. 

4.2. Comparison with previous work 

The simulations presented here suggest that the different dy- 
namics within the plunging region, obtained by Hawley & Kro- 
lik (2001) and Armitage, Reynolds & Chiang (2001), are plau- 
sibly due to different values in the saturation levels of magnetic 



fields in the disk. The relatively strong magnetic fields in the 
former simulations, which obtained f3 5 — 10, are more com- 
parable to the boosted simulations in the current work, and lead 
to substantial angular momentum transport within the plung- 
ing region. Weaker disk fields, with /? ^ 20 and a Shakura- 
Sunyaev (1973) a parameter of a ~ 0.02 — 0.04, produce much 
less striking effects. Although there are undoubtably other sig- 
nificant differences between the simulations, all work to date 
is broadly consistent with the hypothesis that significant mag- 
netic coupling to the plunging region requires disk fields with 
a ~ 0.1 (and /3 ^ 10). This appUes for disks that are mod- 
erately thin, with relative thickness h/r 0.1, which are the 
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only systems that have been simulated so far. 

4.3. Scaling to thinner disks 

Very thin disks are difficult to simulate, and it is not obvi- 
ous how the threshold f3 seen in the current work scales with 
h/r. We note, however, that the Alfven speed in the disk scales 
as, VA oc (^) /3~^/^. Unless /3 in the disk decreases rapidly 
with decreasing h/r, the Alfven speed in the gas immediately 
outside the plunging region will decrease with decreasing disk 
thickness. Conversely, the inflow speed Vr within the plung- 
ing region is a fixed function of r, and is independent of h/r. 
Hence, a larger increase in the Alfven speed is needed if thin- 
ner disks are to maintain causal contact into the plunging re- 
gion. It is therefore probable that it becomes increasingly more 
difficult to extract energy from the plunging region as the disk 
becomes thinner, although further simulation work is required 
before drawing robust conclusions. 

5. CONCLUSIONS 

In this letter, we have used MHD simulations to investigate 
the effect of the magnetic field strength in the body of the ac- 
cretion disk on the dynamics of the material within the plung- 
ing region. For a disk with a sound speed corresponding to 
an h/r « 0.1, we find that plausible saturation values of the 
magnetic field are within a factor of a few of a threshold be- 
tween two regimes. In the runs which saturated with a fairly 
low magnetic field (j3 ^ 20), there was only modest extrac- 
tion of angular momentum from material within the plunging 
region, and the ZTBC was a reasonable approximation. In 
our 'boosted' runs, where the field saturated at a higher level 
(/3 ~ 5 — 10), we found significantly stronger angular mo- 
mentum transport within the plunging region. Furthermore, our 



boosted run showed sporadic intervals of very efficient angular 
momentum transport (with stresses that remained fairly con- 
stant well within the plunging region). We note that stratified 
simulations by other authors (Hawley 2000; Miller & Stone 
2000; Hawley & Krolik 2001) suggest that magnetic fields 
could be stronger (i.e. lower (3) above the disk midplane. Such 
fields, which are not included in our simulations, could con- 
tribute to more significant transport within the plunging region. 

The extraction of energy and angular momentum from ma- 
terial in the plunging region can significantly increase the effi- 
ciency of accretion above the standard value, especially in the 
innermost regions. This raises the possibility of varying effi- 
ciencies, both across systems, and within a given system at dif- 
ferent times. 

In this paper, we have used a large scale flux threading the 
accretion disk as a numerical device for altering the satura- 
tion level of the MRI. Large scale fields could also, of course, 
be present in real systems, and are often invoked in models 
for the formation of jets and outflows (e.g. Ouyed, Pudritz & 
Stone 1997, and references therein). This raises the interesting 
possibility that periods of increased luminosity (due to a high 
radiative efficiency) may be related to periods during which 
strong outflows occur. Such associations have been reported 
for Galactic Black Hole Candidates such as GRS 1915H-105 
(Mirabel et al. 1998). 
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